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We present a fully self-consistent combined GW and dynamical mean field (GW+DMFT) study 
of the spectral properties of the extended two-dimensional Hubbard model. The inclusion of the 
local dynamical vertex stemming from the DMFT self-energy and polarization is shown to cure 
the problems of self-consistent GW in the description of spectral properties. We calculate the 
momentum-resolved spectral functions, the two-particle polarization and electron loss spectra, and 
show that the inclusion of GW in extended DMFT leads to a narrowing of the quasi-particle width 
and more pronounced Hubbard bands in the metallic regime as one approaches the charge-ordering 
transition. Finally, the momentum-dependence introduced by GW into the extended DMFT de- 
scription of collective modes is found to affect their shape, giving rise to dispersive plasmon-like 
long-wavelength and stripe modes. 



PACS numbers: 71.27.+a,72.15.Qm,75.20.Hr 

Modern spectroscopic techniques are able to mea- 
sure one- and two-particle spectra of condensed mat- 
ter systems with remarkable precision, characterizing not 
only quasi-particle excitations but unveiling also satellite 
structures. Examples include Hubbard bands in pho- 
toemission spectroscopy, stemming from the atomic-like 
behavior of the electrons in partially filled narrow d- or 
/-shells p], or collective excitations such as plasmonic 
features. Addressing such effects requires an accurate 
description of one- and two-particle spectral functions 
within the framework of many-body theory. The quanti- 
tative prediction of satellite features has even been used 
as a quality marker for many-body techniques. The fail- 
ure of self-consistent perturbation theory in the screened 
Coulomb interaction, the self-consistent GW approxima- 
tion, to describe plasmon satellites in the electron gas for 
example has provided arguments in favor of a non-self- 
consistent ( "one-shot GW" ) treatment [21 [3] . For real 
solids, few fully self-consistent calculations are available 
[H [5], and no consensus concerning the virtues of self- 
consistency has been reached so far. A popular scheme, 
dubbed quasi-particle self-consistent GW (QPSC-GW) 
[6] yields reasonable estimates both for total energies and 
spectra. Yet, most of the calculations within this scheme 
were applied to semiconductors, and applications to cor- 
related metals only start to appear [7]. The inclusion 
of an appropriate vertex correction is expected to resolve 
the ambiguities around the self-consistency question, and 
it has been in particular proposed that a combined GW 
and dynamical mean field scheme [5] would enable self- 
consistent calculations even for spectral properties. Early 
pioneering calculations on a three-dimensional extended 
Hubbard model EI] have benchmarked several flavors 



of combined schemes along these lines. However, the nu- 
merical difficulty of solving the DMFT equations with 
frequency-dependent interactions has until recently pre- 
vented the direct investigation of spectral properties. 

Implementing the combined GW+DMFT scheme in a 
fully self-consistent manner for the two-dimensional ex- 
tended Hubbard model, we here demonstrate that this 
technique indeed successfully overcomes the deficiencies 
of GW. The implicit inclusion of a non-perturbative lo- 
cal vertex enables fully self-consistent calculations for 
spectral properties. In the correlated metal regime, the 
GW+DMFT self-energy encodes both, band renormal- 
ization effects and Hubbard satellite features. As ex- 
pected from the physical ingredients, the theory also de- 
scribes the Mott insulating state for strong local Coulomb 
interaction, which is inaccessible in GW alone, as well 
as the charge-ordered state driven by intersite interac- 
tions, absent from standard DMFT. In addition, dynam- 
ical screening effects give rise to plasmonic features in 
the local spectral function. While the local spectral func- 
tions in the intermediate to strong correlation regime are 
little affected by the non-local self-energy contributions 
stemming from the GW approximation, a substantial k- 
dependent modulation of the peak widths is observed 
in the momentum- resolved spectral functions. We an- 
alyze momentum-resolved two-particle spectra and show 
that the self-consistent combination of GW and EDMFT 
strongly affects the shape of collective modes, giving rise 
to dispersive plasmon-like long-wavelength modes and 
stripe modes. 

We consider the half-filled extended Hubbard model 



2 



on a two-dimensional square lattice 



H = -t 4a c i<y + J2{Un it n ti - ^n,) + — ^ mrtj, 



where Ci a and c\ a denote the annihilation and creation 
operators of a particle of spin a =\, l at the lattice site i, 
n la = cl a c ia , and n* = n^ + n^. is tne sum over a11 

nearest-neighbor sites, t > is the hopping amplitude be- 
tween two neighboring sites, [i is the chemical potential, 
U the on-site repulsion between electrons of opposite spin 
and V the repulsion between two electrons on neighbor- 
ing sites. The Fourier-transformed bare interaction term 
thus reads = U + 2V (cos(k x ) + cos(fc y )). All energies 
are given in units of the half-bandwidth D = At. We show 
results for inverse temperature /3D = 100, restricting our 
study to the paramagnetic phase. 

The GW+DMFT approach is derivable from a free en- 
ergy functional [1 1 . The Legendre transform of the free 
energy with respect to the Green's function G and the 
screened interaction W can be expressed as a sum of the 
Hartree-Fock part and a Luttinger- Ward-like correlation 
functional ^[G, W], which sums up all skeleton diagrams 
built from G and W [12] • The GW+DMFT scheme con- 
sists in approximating * as * « vE ,EDMFT [G l i, Wu] + 
^nanioc[Gij>Wij], where the first term is calculated from 
a (dynamical) impurity problem as in extended dynami- 
cal mean field theory (EDMFT) [T3HT5] and the second 
term is the non-local (i ^ j) part of the GW functional 

*^oc[G«, Wij\ = E w (■:.;".,(■:„■ 

The GW+DMFT scheme self-consistently constructs 
the Green's function G and the screened interaction 
W of the system as a stationary point of the free en- 
ergy functional. The self-energy £ and polarization P 
are formally obtained by functional differentiation of 
with respect to G and W, leading to the expressions 
E(k,ioj) = Y, imp (iuj n ) + S°^ DC (fc, iu) and P{k,iv n ) = 
Pimp(w n ) + Pn^\oc( k ^ iv n) ( w n an d v n are fermionic and 
bosonic Matsubara frequencies, respectively). This en- 
dows GW+DMFT with conserving properties [T6 . The 
momentum-dependent G and W are then calculated from 
the one- and two-particle Dyson equations and used as 
inputs for a GW calculation, yielding £ GW and P GW . 
Meanwhile, their local parts are extracted to compute 
the local Weiss fields Q and U: Q^ x (iu n ) — G^ c {iu> n ) + 
S imp (iw„) -atlAU^ 1 (iv n ) = W^ l(w n )+P- imp (iv n ). These, 
in turn, are used as inputs to a dynamical impurity 
model, which we solve using a continuous-time Monte 
Carlo algorithm [TTJ [Tg] to obtain updated local self- 
energies. The whole scheme is iterated until convergence. 
The calculations have been performed on a 64 x 64 mo- 
mentum grid, while the analytical continuation of the 
imaginary-time data has been performed using the Max- 
imum Entropy method [TS]. We monitor the follow- 
ing quantities: (i) the local spectral function A\ oc {u) = 
— -ImGioc(w), (ii) the local and non-local self-energy, 
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Figure 1: Spectral function Ai oc (uj) obtained within four 
different self-consistent schemes (see text). 



(iii) the local and non-local polarization, (iv) the electron 
energy-loss spectrum (EELS) Im [— e(k, w) _1 ] (where e is 
the dielectric function, e(k, lj) = 1 — VkP(k, oj)). 

Within extended DMFT and GW+DMFT, in the ab- 
sence of intersite repulsion, the Mott transition takes 
place at U c « 2.5. This value is slightly modified by 
intersite repulsions V < V c = 0.8. At V c a transition to a 
charge-ordered phase occurs. In the following we study 
the local spectral properties in the metallic phase with 
weak (U = 0.5, V = 0.1) and intermediate (U = 1.5- 
2, V = 0.4) interactions as well as in the Mott insu- 
lator at U — 3.5 and V = 3. Figure [T] shows the lo- 
cal spectral function A\ oc (uj) obtained within (i) (self- 
consistent) EDMFT, (ii) self-consistent GW+DMFT, 
(iii) self-consistent GW and (iv) quasi-particle self- 
consistent GW (QPSC-GW). The latter scheme was im- 
plemented by computing the lattice Green's function 
from the GW self-energy via G(k, iw„) _1 = iu n — Zk{ek — 
Re£ GW (fc,Kj )), where Z k « (1 - ImS GW (fc, iw )/wo) _1 
is the quasi-particle weight as estimated from the value 
of the self-energy at the first Matsubara frequency. For 
small interactions (panel la), correlation effects are neg- 
ligible, and the four schemes result in indistinguishable 
spectra within the numerical accuracy. As the local inter- 
action U becomes significant (panel lb) deviations start 
to appear, with EDMFT and GW+DMFT exhibiting 
stronger correlation effects. Upon increasing local in- 
teractions (panel lc), the quasi-particle renormalization 
becomes stronger, the width of the coherent central peak 
shrinks, and the corresponding spectral weight is trans- 
ferred to higher energies. This - physically expected - be- 
havior is realized by the EDMFT and GW+DMFT spec- 
tra, which exhibit higher-energy structures at oj rs +C//2 
already for U = 1.5. The Hubbard bands gain spec- 
tral weight as U increases further. Interestingly, they 
are more pronounced within GW+DMFT than within 
EDMFT. Finally, at U = 3.5, a Mott gap has opened, 



3 





1.2 



EDMFT 
GW+DMliT 



1.5 



1 



0.8 - 



A(x,ji)(«>) 
0.4 U= 2 . V = 2 



EDMFT 
GW+DMFT 

0.8 - 
0.5 A M (co) 
0.4 -U=2, V=.4 




(0,0) (Ji,7i) (71,0) (0,0) 



Figure 2: Spectral function Ai oc (u;) obtained using 3 different 
one-shot GW schemes (see text). Lower right panel: vertex 
estimate A(uj = 0) as a function of U. 



and EDMFT and GW+DMFT spectra are similar. In 
addition to the two Hubbard bands, the EDMFT and 
GW+DMFT spectra display two symmetric high-energy 
satellites, whose spectral weight depends on the intersite 
interaction V. The QPSC-GW spectra display only a 
weak renormalization of the bandwidth as U increases 
from the weak to the strong coupling limit, and at 
all correlation levels the spectra remain metallic. The 
same is true within the self-consistent GW method, al- 
though with increasing correlations some spectral weight 
is shifted to higher frequencies, albeit in a featureless way. 

These observations show that both self-consistent GW 
approaches yield a correct result only in the weak- 
correlation regime. As correlations increase, GW fails to 
describe the shift of spectral weight to high-energy inco- 
herent bands, present in DMFT. We note that in the local 
GW+DMFT spectra the Hubbard bands are enhanced 
compared to the EDMFT or GW spectra. This effect 
can be ascribed to the self-consistency, which allows the 
local quantities to re-adjust to the non-local self-energies 
E G w and P G w- 

Another salient characteristic of EDMFT and 
GW+DMFT spectra is the presence of additional high- 
energy satellites in the Mott phase. These directly reflect 
the frequency-dependence of the local interactions U (cj) 
induced by the nearest-neighbor repulsion term V . As 
demonstrated in Ref. [2U], a pole in U{d) (such as a plas- 
mon pole) leads to multiple satellites in the local spec- 
tral function. In our case, the local interactions U(oj) are 
characterized by a broad continuum of poles centered at 
some energy Wd, resulting in only two symmetric satellites 
in the Mott spectra. In the metallic phase, these satel- 
lites are present, but they are broad and merged with the 
Hubbard bands, making them hardly distinguishable. 

The failure of both self-consistent GW schemes to cap- 
ture Hubbard bands or high-energy satellites is consistent 



Figure 3: Left panels: momentum-resolved spectral func- 
tion A(k, u) at U = 2, V = 0.4 within EDMFT (top) and 
GW+DMFT (bottom). Right panels: A(k,u) at k = (ty,tt) 
for U = 2, V = 0.2 (top) and U = 2, V = 0.4 (bottom) within 
both schemes. 



with the well-known observation that self-consistency in 
GW for the homogeneous electron gas results in a smear- 
ing out and displacement of high energy satellite fea- 
tures [2]. In light of this observation, most modern 
GW schemes therefore adopt a "best-G-best-W" strat- 
egy, rather than aiming at full self-consistency. Figure 
2 illustrates the virtues and limitations of this strategy 
by displaying the spectra obtained in different one-shot 
GW schemes: (i) in "GqW" , the non-interacting Green's 
function G and the converged GW+DMFT W are taken 
as inputs to a one-shot GW calculation, (ii) "GWo" takes 
the converged GW+DMFT G and W = v(I-vGqGq)- 1 
within the random phase approximation as inputs and 
(iii) "best G, best W" takes the converged GW+DMFT 
G and W as inputs. At all correlation levels (U = 0.5 
to U — 3.5), these three GW schemes produce results 
very similar to self-consistent GW. In particular, they 
remain metallic. Even the "best G, best W" scheme 
in the Mott phase (U = 3.5) yields a metallic self- 
energy, despite the Mott-like character of the input G 
and W. This phenomenon is due to the lack of Hedin's 
three-legged vertex A in GW schemes, as shown in the 
lower-right panel of Figure 2. There, an estimate of 
the local part of A is computed from EDMFT results 
at V = 0. Remembering that, schematically, the irre- 
ducible vertex function A appears in the self-energy as 
£ = GWA [5T], a rough estimate - neglecting the true 
frequency structure - is computed as follows: one com- 
putes an effectively vertex-corrected screened interaction 
W(t) = S lmp (r)/G imp (r) from EDMFT, then Fourier- 
transforms it to W(iv n ); finally, the static vertex esti- 
mate is obtained as A(0) « W(«i / o)/Wioc(« l 'o)- Crude as 
it is (the full vertex depends on two independent frequen- 
cies), this estimate nonetheless clearly demonstrates the 
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Figure 4: U = 2, V = 0.4. Upper panels: ImP(fc,w) 
within EDMFT (left) and GW+DMFT (right). Lower panels: 
Ime _1 (A:,w) within EDMFT (left) and GW+DMFT (right). 

role of vertex corrections for the Mott transition: from 
unity in the weakly correlated regime, it increases with 
U until it diverges at the Mott transition. This indicates 
that within the language of Hedin's equations, the diver- 
gence of the local vertex is the driving force of the Mott 
phenomenon, making any vertex-less approximation un- 
fit to capture it. 

The effect of the non-local GW contributions on 
EDMFT are illustrated by the momentum-resolved spec- 
tral functions, displayed in Figure [3] Generally, the dis- 
persion of the Hubbard bands follows the dispersion of 
the quasi-particle peak, within both schemes. In the pres- 
ence of a strong intersite interaction, the non-local self- 
energy and polarization lead however to an additional 
/c-dependent modulation of the linewidth and weights. 
This is illustrated in the lower panels of Fig. [3] for the 
(7r, 7r) point, where a pronounced sharpening of the quasi- 
particle peak is observed along with enhanced weight of 
the Hubbard bands. 

We now turn to a study of two-particle quantities. 
Figure [4] shows the momentum-resolved imaginary part 
of the polarization and the electron energy loss (EELS) 
spectrum Im [— e _1 (fc, w)] in the metallic regime. Within 
EDMFT, the polarization displays a broad mode which 
reflects the particle-hole excitations of the system. They 
are centered at U/2, reflecting the emergence of the 
Hubbard bands and the corresponding excitations be- 
tween Hubbard bands and the quasi-particle peak. In 
contrast, the polarization spectrum within GW+DMFT 
is dispersive. While displaying sharper features close 
to the r = (0, 0) point, it captures particle-hole ex- 
citations due to Fermi-surface nesting at wave-vector 
(7r,7r), as well as the zero-sound mode at long wave- 
lengths and low energies. The EELS spectrum contains 
the particle-hole excitations (poles of the polarization) of 
the system and its collective modes, which correspond 



to the solutions of TLeP(k,ui) — \/vt- These collec- 
tive modes are damped out close to particle-hole exci- 
tations (when ImP(fc,w) is large). This analogue of the 
free-electron-gas Landau damping occurs at the (7r,7r) 
point in EDMFT and GW+DMFT. It can be directly as- 
cribed to the nearest-neighbor repulsion, which induces 
scattering along this direction. The energy and life- 
time of this collective excitation differs from EDMFT to 
GW+DMFT. In GW+DMFT it is lower in energy, more 
dispersive and with a larger lifetime. In GW+DMFT, 
two modes are visible above the (7r,0) point, indicat- 
ing the existence of two stripe modes at energies oj = 1 
and u> — 2 corresponding to stripe-like modulations, 
where the sign of the density fluctuation varies from row 
to row in the x-direction. For obvious reasons, they 
are not captured by EDMFT. These two-particle excita- 
tions are directly related to the screening in the system 
as the screened interaction, W, is given by W(k,oj) — 
e~ 1 (fc,w)ufc. In particular, they explain the retardation 
effects at play in the local interactions U(uj) and the cor- 
responding satellites in the local spectra. 

In conclusion, we have demonstrated how the am- 
biguities on the optimal degree of self-consistency in 
many-body perturbation theory are resolved by includ- 
ing a non-perturbative local vertex in the calculation. 
Based on a fully self-consistent implementation of the 
combined GW+DMFT scheme, we have analyzed one- 
and two-particle satellite features in correlated materi- 
als. While we confirm the well-known "washing out" 
of satellite features in self-consistent GW calculations, 
self-consistent GW+DMFT docs not suffer from this de- 
ficiency. Plasma- and zero-sound-like oscillations involv- 
ing itinerant carriers as in the electron gas survive only 
in the regime of small local Coulomb interactions, but 
are quickly suppressed in the correlated metal. In this 
regime, excitations related to the creation of doublons 
become dominant. While local spectral functions are lit- 
tle affected by non-local contributions in a wide range 
of parameters, the momentum-resolved spectra display a 
fc-dependent modulation of the width of the peaks; the 
momentum-dependence introduced by the GW part be- 
comes truly crucial when assessing dispersions of two- 
particle spectral properties, differentiating in particular 
the nature of the collective modes in the (0,0), (0,7t), 
and (tt,it) directions. Our findings have implications for 
the nature of satellite features in correlated materials. 
In particular, it becomes obvious that electron-gas-like 
plasmons in materials stem dominantly from the charge 
contained in completely filled shells (that is from multi- 
orbital effects), while partially filled shells give rise to 
doublon excitations of the kind we describe. The inter- 
play of local correlations and charge ordering phenomena 
and their intriguing wave-vector dependence call for an 
extension of our study to realistic two-dimensional sys- 
tems: recent experimental findings of charge ordering in 
cuprates [35] , cobaltates [23J or in systems of adatoms on 
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surfaces [Uj are prominent examples. 
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